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ABSTRACT 

Context. The existence of rapidly rotating cool stars in young clusters implies a reduction of angular momentum loss rate for a certain 
period of the star's early life. Recently, the concentration of magnetic flux near the poles of these stars has been proposed as an 
alternative mechanism to dynamo saturation in order to explain the saturation of angular momentum loss. 

Aims. In this work we study the effect of magnetic surface flux distribution on the coronal field topology and angular momentum loss 
rate. We investigate if magnetic flux concentration towards the pole is a reasonable alternative to dynamo saturation. 
Methods. We construct a ID wind model and also apply a 2-D self-similar analytical model, to evaluate how the surface field 
distribution affects the angular momentum loss of the rotating star. 

Results. From the ID model we find that, in a magnetically dominated low corona, the concentrated polar surface field rapidly 
expands to regions of low magnetic pressure resulting in a coronal field with small latitudinal variation. We also find that the angular 
momentum loss rate due to a uniform field or a concentrated field with equal total magnetic flux is very similar From the 2D wind 
model we show that there are several relevant factors to take into account when studying the angular momentum loss from a star 
In particular, we show that the inclusion of force balance across the field in a wind model is fundamental if realistic conclusions are 
to be drawn from the effect of non-uniform surface field distribution on magnetic braking. This model predicts that a magnetic field 
concentrated at high latitudes leads to larger Alfven radii and larger braking rates than a smoother field distribution. 
Conclusions. From the results obtained, we argue that the magnetic surface field distribution towards the pole does not directly limit 
the braking efficiency of the wind. 

Key words, stars: late-type - stars: magnetic fields - stars: winds, starspots 



1. Introduction ena occurring in the stellar corona. It can cause non-solar 
5— i type phenomena like flares and X-ray emission at high lat- 
d Rapidly rotating cool stars have a surface magnetic flux dis- jtudes jS chmitt & Favata, 1999) and sUngshot prominences 
tribution significandy diff^erent from that observed on the Sun. dCoUier Cameron & Robinsoft,1989D . Here, we concentrate our 
Doppler imaging, and more recently, Zeeman Doppler imag- attention on its imphcations for steUar spindown. 
ing, gives clear evidence for spot s and signific ant surface mag- 
netic flux at high latitudes (e.g. | Donati etaU 1999). In many ^ind magnetic braking is based on the principle that when 
cases, a large spot, or cluster of spots^s present at the pole emitted from a star is kept co-rotating with die star by 
d Vogt & Penrodj [19831 I Strassmeieij |2002i). But clearly, spots magnetic torque, it transports significantly more angular mo- 



and magnetic fields are also present at low latitudes (e.g. mentum outwards than gas that conserves its angular momen 
[ Barnes etaU 1998). Knowledge of the large-scale coronal mag- as it moves outwards (ISchatzmanL \l96^ . In axisymmet- 
netic topology can in principle be obtained b y extrapolating winds the angular momentum loss rate is equivalent to that 
magnetic fields from the boundary data (e.g. [Hussain et al.|, carried by the gas kept in co-rotation with the star out to the 
2001). However, the absence of information from a fraction of the Alfven surface, where the poloidal wind velocity equals the 
the stellar surface and the weak coiTelation between spots and poloidal Alfven velocity. The presence of rapidly rotating late- 
magnetic signatures makes us look on the results of this promts- jypg gj^rs in young clusters implies that there must be some 
ing technique with some caution. limitation to the eflicienc y of magnetic br aking during the pre- 
The presence of spots and surface magnetic flux at high main sequence phase (e.g. lBarnes & So"fial 1996). One possibil- 
latitudes has important consequences for several phenom- ity is that beyond some rotation rate the magnetic field strength 

no longer incre ases with increasing rotation r ate, i.e., the dy- 

Send offprint requests to: A. Aibeo, e-mail: namo saturates dMacGregor & Brenneil Il99lh . But several al- 

aaibeo@demgi . estv. ipv.pt ternatives to dynamo saturation have been suggested. These in- 
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elude, the increase of t he closed field region with rotation rate 
dMestel & Spruid 1 1987b . and the concentration of magnetic flux 
near the pole s of rapidly rotating stars (Solanki et al., 1997; 
iBuzasiL fT997h can lead to a saturation in the angular momen- 
tum loss rate. The effect of the closed field region, or dead zone, 
on the rate of braking has been stud ied in some detail for both 
single and binary stars (e.g. lLi et all 1994). 

The idea that the concentration of magnetic flux near the 
poles can mimic dynamo saturation is based on basic princi- 
ples. The Alfven radius is large near the poles where the field is 
strong, but smaller near the equator where the field is weaker As 
wind braking is mainly from the contributions of low and inter- 
mediate latitudes, the overall effect is a reduction in the braking 
when compared with an homogeneous surface field distribution 
wit h the same total magn etic flux. In recent studies, the classi- 
cal l^^bgE&D^ii (Il967h wind model is explored in detail and 
extended to study the effect of non-uniform surfa ce magnetic 
field distribution on the angular m omentum loss rate jHolzwarthl 
l2005l:lHolzwarfli & Jardin^l2005h . All these models assume that 
the coronal radial field has a latitudinal distribution similar to the 
surface radial field, i.e., the force balance across the field is ne- 
glected. 

In the particular case of the Sun, Ulysses observations have 
shown that there is no significant gradient in latitude in the radial 
component of the interplanetary magnetic field (Balog h et al.L 
ll995l:ISmith & Balogh, 1995). These observations have been ex- 
plained as a result of the low plasma beta of the solar corona 
so that the latitudinal and longitudinal gradients in radial fields 
rel ax quickly, creating a n essentially uniform field by 5-10 
R^ (ISuess&SmithLlT996h . 

In the present work we consider whether polar magnetic flux 
concentration towards the pole is a valid alternative to dynamo 
saturation as an explanation for angular momentum loss satura- 
tion. The aim of this paper is also to demonstrate that the inclu- 
sion of force balance across the field in a wind model is funda- 
mental if realistic conclusions are to be drawn from the effect 
of non-uniform surface field distribution on magnetic braking. 
Section|2]presents a ID Weber & Davis type wind model along 
a totally opened potential magnetic field resulting from a surface 
flux distributio n concentrated tow ards the pole. In Sect.[3]the 2D 
wind model of iLima et al.l (120011) is used to study the effect of 
the variation of surface flux with latitude on the rate of magnetic 
braking. The implications and limitations of the results obtained 
are discussed in Sect.|4]and the conclusions presented in Sect. El 



2. Application of a 1 D wind model 

Present models of the inhibition of the angular momentum loss 
by polar concentration of the magnetic field are quahtative in 
nature and ne glect fo rce balance across the field (ISolanki et al.L 
119971: iBuzasil 119971 iHolzwarflj, l2005h . They are based on the 
assumption that the concentration of magnetic field at high lati- 
tudes generates a large Alfven radius near the pole and a small 
Alfven radius near the equator leading to a smaller effective 
Alfven radius and therefore a reduced braking efficiency. 

Ignoring force balance across the field and the indirect ef- 
fect of the magnetic field on the wind dynamics are too severe 
restrictions present in these qualitative models. This makes us 
doubt the conclusion that field concentration at high latitudes 
significantly reduces the effective Alfven radius and angular mo- 
mentum loss rate. Therefore, w e present a simple wind model 
akin to these qualitative models dWeber & Davisi. [19671: ISakurail 
119851: iHolzwarflil l2005h for a surface field concentrated towards 



high latitudes but where force balance across field lines is par- 
tially taken into account. In modelling the solar corona and so- 
lar wind it is common to apply the po tential field source sur- 
face model of the coronal ma gnetic field jAltschuler & NewkirkL 
EM Is uess & SmifliL[T996l) . In this model, the magnetic field is 
assumed to be potential between the surface and a spherical outer 
surface where the field is required to become radial. Although 
this model ignores volume and surface currents, it would be ad- 
equate to show that the latitudinal gradients of the radial field 
are s moothed ou t at relativ ely shor t distances from the surface 
(e.g. lSuess et all 1977; Rile v et aTl 2006). However, this model 
would be inadequate to determine the field configuration outside 
the source surface and useless if one assumes the field to be to- 
tally open. Therefore, we use a different model to determine the 
coronal magnetic field. 

We construct a totally open magnetic field configuration by 
considering the magnetic field to be the dominant force in the 
low corona so that it is potential everywhere except at an equa- 
torial current sheet. We then solve the wind equations along the 
field lines and compare our solutions with those that result from 
a uniform surface field. 

We start by considering an axisymmetric poloidal magnetic 
field in an atmosphere with negligible mass and pressure. To al- 
low for the effect of the wind without directly solving the com- 
plicated set of equations, a certain amount of magnetic flux is 
taken as open. We then make use of a family of analytical so- 
lutions to construct axisymmetric, partially open magnetic fields 
that are pote ntial everywhere except on a force -free equatorial 
cun-ent sheet (! Lowlll986l;lLepeltier & Alvlll996l) . The magnetic 
field can be expressed in spherical coordinates in terms of the 
stream function A: 
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so that Maxwell's equation V.B = is satisfied. In this way, 
magnetic field lines are represented by contours of constant val- 
ues of the stream function, A. The stream function is given by a 
linear combination of basic functions Z„ 



(2) 



where n takes odd values, y„ are constant coefficients and the 
analytical functions Z„ are developed from the oblate spheroidal 
harmonics and classified according to the harmonic order n. 
There is a free parameter, a, representing the radial distance 
beyond which the magnetic field is completely open. Because 
these functions are not orthogonal, it is complicated to deter- 
mine the coefficients y„ for a prescribed boundary condition. 
Nevertheless, we can combine different functions with suitable 
coefficients to obtain a field with the desired properties. The pro- 
cedure to determine the different Z„ is described in Low (1986) 
where Z\ an Z3 are explicitly given. The functions Z5 and Z7 are 
given in AppendixlAl We consider the case Adip - yf ''Zi , repre- 
senting a dipole-like field and A^o/a,. = yiZ\ +yT,Z^ +JiZ=, +y-iZ-i, 
representing a magnetic field concentrated at the poles. The con- 
centrated field distribution is not intended to be a realistic field 
distribution but merely an extreme case of field concentration 
towards high latitudes used for illustrative purposes. 

We consider fully open fields in accordance with the pic- 
ture presented in previous works (^Solanki et al., 1997; Buzas^, 
Il997i iHolzwarth & JardinH |2005|) . Thus, we set a = ro (where 
ro represents the stellar radius) and use the values of y„ given in 
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Table 1. Coefficients y,, for a fully open magnetic field 4 



Fully Open 

n 1.0 

73 -3.2 X 10-2 
75 1.4 xlQ-^ 
77 1.0 xlO-'* 
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Fig. 1. Fully open fields. The radial field latitudinal distribution: 
at the surface (full line), at r = A-ro (short dashed) and at r = 8ro 
(long dashed). The thick lines represent the poleward concen- 
trated field and the thin lines the dipole-like field. For a clearer 
representation all the radial fields are set to unity at the pole. 



Table [T] while the value of yj''' is determined by the condition 
that the total magnetic flux is the same in both cases. 

Figure [1] represents the radial field strength as a function of 
co-latitude for different radial distances from the stellar center 
for the two cases considered. The latitudinal profiles of the two 
radial fields are very different from each other near the surface 
of the star. However, at intermediate (r=4ro) and at large dis- 
tances from the surface (r=8ro) they are similar to each other 
and to the split monopolar field. In Fig.|2]the coronal field topol- 
ogy for both cases is represented. This figure clearly shows that 
for surface fields concentrated near the poles, the field lines that 
emerge at high latitudes are pushed towards low latitudes in the 
low corona, resulting in a field almost independent of latitude 
further out in the corona. It is straightforward to show, using the 
same model for a > ro, that the same happens for partially open 
fields. 

These results demonstrate that very different surface flux dis- 
tributions give rise to similar coronal fields. Physically, this is 
simply a consequence of the very large magnetic pressure dif- 
ference between high and low latitudes (this is clearly illustrated 
in Parker's monograph, Parker ( 1979)). Therefore, a strong field 
near the pole and a weak field near the equator at the stellar sur- 
face do not imply that the same is true in the corona. We note 
that, in essence, this argument is identical to the one used to 
explain Ulysses observations concerning the lack of latitudinal 
gradi ents in the radial component of the solar wind magnetic 
field dSuess & Smithlll996i) . Having obtained the poloidal field 
configuration, we now determine the polytropic wind solution 
along the field. We consider a star of one solar mass and ra- 
dius, rotating rigidly and characterized by a corona with uni- 
form base temperature, T - 2.1 x lO^K, and uniform base den- 
sity, p = 10"'^ kg/m^. We also assume a mean atomic weight 
of yu = 0.6 and a polytropic index of 7 = 1.15. Here we fol- 
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Fig. 2. Fully open fields. The lines of force of the dipole-like 
surface flux distribution (dashed line) and the flux distribution 
concentrated towards the pole (full line). For the later case notice 
the field lines bending towards the equator. 



low closely the approach described in ISakurail (Il985 [). In brief, 
we solve the Bernoulli equation presented in the Appendix |B] 
for different stellar rotation periods, ranging from 1 to 30 days, 
and for different total surface magnetic fluxes, corresponding to 
uniform fields ranging from 1 to 100 G. However, we do not 
assume the area of each flux tube to increase as r^, but we im- 
pose the area variation to be that of the solution of the field. In 
particular, we consider two flux tubes with contrasting proper- 
ties. One has its roots near the pole at latitude 89° where the 
field near the surface expands very rapidly, and the other has its 
roots close to the equator at latitude 1 ° where the field near the 
surface first contracts and then expands very slowly. Both flux 
tubes expand oc far from the surface. The latitude of the flux 
tube is assumed the same and constant for all the flux tubes. We 
then compare the solutions obtained with those along a flux tube 
with an area variation oc r^. The solutions are compared impos- 
ing equal flux tube area at the Alfven radius as the angular mo- 
mentum loss rate is determined from the Alfven radius, velocity 
and density (Eq.fTSI). Density, velocity and radius at the Alfvenic 
point are denoted by p», v» and respectively. We find that: i) 
the Alfven radius changes very little (< 2%), being larger for 
rapidly expanding loops; ii) the base velocity increases consid- 
erably as one goes from slowly expanding to rapidly expanding 
flux tubes (over two order of magnitude); iii) the mass loss rate 
at the Alfen radius changes little from flux tube to flux tube; iv) 
the slow magnetosonic point decreases significantly from slowly 
expanding to rapidly expanding loops (up to 25%); v) the differ- 
ences in angular momentum carried by the wind and magnetic 
torques (oc p^VtR*) by the different flux tubes of equal area at R, 
are negligible (< 2%). These results determine that, in disagree- 
ment with previous results, there is no significant difference in 
the angular momentum loss when the surface field changes from 
being uniform to being concentrated towards the poles. We note 
that centrifugal acceleration on rapidly rotating stars is ineffi- 
cient at high latitudes and generates a non-spherical Alfven sur- 
face, but clearly, this is independent of the surface field distribu- 
tion. 

In the simple model presented here we observe the expan- 
sion of the polar field to low latitudes. However, this expansion 
can be limited by a high equatorial gas pressure and the Lorentz 
force due to the azimuthal magnetic field. Therefore, these par- 
ticular results are only valid provided the corona near the stel- 
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lar surface is a low beta plasma and the effect of the azimuthal 
magnetic field can be neglected, which is likely to be the case 
in magnetized rapidly-rotating stars. The parameters considered 
range from those characteristic of stars termed as slow magnetic 
rotators (SMR) to thos e characteristic of fast magne tic rotators 
(FMR), as classified bv lBelcher & MacGregoJ(ll976l) . The more 
rapidly rotating and magnetic stars considered in the present 
analysis are F MR, for which there is a pole ward collimation of 
the field lines dHevvaerts & Normani Il989h that is neglected in 
these simple models. But, as will be argued in detail in Sect. |4] 
this is largely independent of what is under scrutiny here - the 
effect of the surface magnetic flux distribution. 

To complete the model just presented, two different and 
complementary approaches could be pursued. In one ap- 
proach, numerical simulations are employed (e.g. ' Sakurail 
1985; [Keppens & Goedbloed, 1999) while in the other, ana- 
lytical or semi-analytical solutions of the equations are ob- 
tained under some simplifying assumptions. Several self sim- 
ilar non-polytropic models have been developed based on 
the assumpti on of a non-li near separation of v ariables (e.g. 
iTsing anos & Trussonil 1991; ISautv & Tsinganosl 1994). Here, 
we adopt the analytical 2-D model of iLima etaP (l200lh as It 
is the only available analytical model that allows us to vary the 
magnetic surface flux distribution and determine how it affects 
the wind dynamics and the angular momentum loss rate. As in 
the 1-D model of Holzwarth (2005), the poloidal field and ve- 
locity are purely radial, but, in this case, force balance across 
the field is obeyed in a self-consistent way. Although these ana- 
lytical models are more rigorous than the qualitative model just 
presented, they require several simplifying assumptions that may 
be physically unrealistic or undesirable. Therefore, it seems nec- 
essary to study how such analytical models compare with quali- 
tative models which are more commonly applied to address the 
problem of magnetic braking. 



3. Application of a 2D wind model 

3.1. A self-similar MHD wind model 

In order to construct a model f or an axisymmetric wind emanat- 
ing from a central rotating star. lLima et al.l (1200 Ih have assumed 
0-self similarity and deduced a solution of the system of ideal 
MHD equations. The solution is found by a non-linear separa- 
tion of variables, keeping the treatment as general as possible 
(e.g. not assuming, a priori, any prescribed variation with lati- 
tude of the velocity, density or magnetic field). 

The model uses spherical coordinates [r, 6, 0] and assumes a 
simple geometry with zero meridional components of the veloc- 
ity and magnetic fields in order to find a treatable form for the 
fundamental solutions an therefore calculate a solution. The out- 
flow dynamics are described by the following set of equations 
of distance and co-latitude for the radial velocity, azimuthal ve- 
locity, radial magnetic field, azimuthal magnetic field, density, 
pressure and stellar angular velocity: 



^Rl/R^ - 1 



B^(R,0) = ABq 



Vr(R,d) = VoY 



1 + n sin^^ e 



+ 5 sin^^ 



V^(R,0) = AVo 



(Y,-Y 



2\ 



R sin" 



Vl +6sm^'0 
B,.{R,0) = ^^1 +fism^'0 



(3) 
(4) 
(5) 



R sin'' 



p(R,0) = ^(l+^sin^^e) 
P(R,0) = ^poV^{Qo + Qism^'0) 



Q(0) = 



AVoY, sin'-^ 



^0 y/l+6sm^'0 



(6) 
(7) 
(8) 
(9) 



where po, Bq, Vq represent the density, radial magnetic field 
and radial velocity at the polar base of the wind, respectively, 
and Ma, Qo and Qi are functions of R and Y, itself a function of 
R. The function Ma is the Alfven number defined by the ratio of 
the poloidal velocity to the Alfven velocity (M^ = AnpV} I B^^). 
From this definition and from Eqs. [3] and |5] it results that the 
Alfven iso-surfaces are spherical. The calculation of the function 
expressing the radial dependence of the radial velocity, Y{R), is 
made from the combination of the radial and latitudinal compo- 
nents of the momentum equation, which results in one first order 
non-linear differential equation. Such an equation, which com- 
bines force balance along and across the fieldlines, shows two 
points where both the numerator and denominator vanish simul- 
taneously. These singular points are related to the non-linearity 
of the steady-state system of equations. The Alfven point, where 
Ma - 1, is a star-type singular point and it is indicated by 
(7?,, F,). All solutions can pass through it. The second singu- 
lar point is the fast magnetosonic point and it is an X-type point 
allowing only two solutions to cross it. Due to the self-similar 
nature of this model, in which force balance is solved simul- 
taneously along and across the fieldlines, the results show that 
there are only the above two critical points. Other self-similar 
wind models shovv the same number of critical points (e.g. 
ISautv & T singanos^. 1994). A detailed discussion of the nature of 
these critical points can be found in Tsinganos et al., (1996). The 
functions Qo and Q\ represent, respectively, the isotropic and lat- 
itudinal dependent normalized compon ents of the p ressure (for 
further details on their calculation see iLima et aL i 2001). The 
temperature is calculated a posteriori using the classical ideal 
gas law. 

The solutions of this model are defined by six dimension- 
less parameters: A, v, M°, ji, 5, e. The parameter A represents 
the ratio between two velocities: the equatorial stellar rotation 
velocity and the polar radial velocity at the Alfven point: A - 
rciQ.{RQ, 7T /2) / VriR*, 0) (for 6 = 0). As a consequence of the an- 
gular momentum conservation and induction equation, this pa- 
rameter can analogously be defined using the two components 
of the magnetic field. The parameter v is the ratio of the stellar 
escape velocity and Vq. Thus, for a given value of v, the param- 
eter A and Y, define how fast or slow rotator the star is. The 
parameter defines the Alfven number at the wind base, i.e., 
it defines how magnetized the star is. The remaining three pa- 
rameters are related to the latitudinal distribution of the different 
physical quantities. 

The parameter 6 evaluates the density anisotropy between 
the equator and the pole (cf Eq.[8]). In a similar way, the param- 
eter fi determines the radial magnetic field anisotropy between 
the equator and the pole (cf. Eq.|5]). Finally, e controls the shape 
of the latitudinal distribution of the magnetic field, density, and 
velocity where high values lead to steep variations and low val- 
ues lead to smoother variations. These three anisotropy parame- 
ters yield the flexibility of generating solutions corresponding to 
different latitudinal dependences and enable us to model stellar 
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R.=17.05 

=260.4 [km s"^] 
R^=24.83 
v'=381.5 [km s"^] 
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Fig. 3. A typical wind solution for a highly magnetized object in the 
iLima et al] l |2001|) model, featuring two distinct singular points; the 
Alfven point and a fast magnetosonic point. 



outflows showing distinct latitudinal distributions of magnetic 
flux. 



3.2. Analysis of a typical wind solution 

From a purely theoretical point of view, if the concentration 
of open magnetic flux reduces the efficiency of magnetic brak- 
ing, then it will do so for both fast and slow rotators. Both ob- 
servations and theoretical modelling indicate that a substantial 
amount of magnetic flux concentrates at high latitu des for solar 
type-s tars of periods below a few day s (~ 5 dav: IStrassmeieii 
2005. ISchuessler & Solankl 1992 and lSchriiver & Titld 2001V 
We then consider a young, moderately rapidly-rotating solar type 
star with ro = 1 R©, M - I Mq, with a 4 day period and with a 
coronal density po = lOp© - 1.6 x 10"'^ kg/m^^. The param- 
eters adopted are: M^" = 0.01, v = 80, /I = 0.11, yielding 
Vo = 7.72 km/s, Vi = 0.85 km/s, and Bq = IIG. The param- 
eters describing the latitudinal profiles are 5 = 5, ju = -0.8 and 
6=1.0 that are representative of a star with more magnetic flux 
at high latitudes than at low latitudes. 

The wind solution obtained is presented in Fig.[3]and in order 
to understand what forces are relevant in accelerating the wind 
and in maintaining the equilibrium in the latitudinal direction, 
we present a detailed study of this solution. 

The three components of the equation of motion, under the 
assumptions of axisymmetry and zero theta component's are ex- 
plicitly written in Eg s . [TOl to [T2l where we have labelled the dif- 
ferent terms with roman numerals for easier identification. 



dVr 



_ dp B^dB^ pGM 

^ r dr Anr An dr 



II III IV 



V 



VI 



jcosfl dp BrdB,. B^cosO B^dB^ 

* sine "dd' A^~de ~ 4^ sinfl ~ J^'dO 
VII VIII IX X XI 



pVr d(rV^) B, d(rB^) 



dr 



Anr dr 



(10) 



(11) 



(12) 



In Fig.|4]we can observe that close to the surface the dom- 
inant forces in the radial direction are the gas pressure gradient 
(III) and gravity (VI), which balance each other in almost hydro- 
static equilibrium. Further out, at high and intermediate latitudes 
(panel a) of Fig.lUi, it is the pressure gradient that accelerates the 
wind with a small contribution from the other forces. However, 
near the equator both the centrifugal force (II) and the Lorentz 
force (IV -I- V) due to the 0-component of the field are important 
in accelerating the wind near and beyond the singular points. 
This is in accordance with what one would expect for a star with 
the properties considered here. 

It is also instructive to analyze the force balance in the 0- 
direction. Figure |5] shows that below the critical points the mag- 
netic pressure gradient due to the radial magnetic field (IX) is 
balanced by the gas pressure gradient (VIII) and, to a lesser ex- 
tent, by the Lorentz force due to the ^-component of the field 
(X-hXI). The gas pressure gradient decreases and changes sign 
close to the Alfven point, as indicated by the crosses on the dot- 
ted curve. In this region the force opposing the radial magnetic 
field pressure gradient is the Lorentz force resulting from the 
0-component of the magnetic field (X + XI). Therefore, the ex- 
pansion of the poloidal field from the high latitudes towards the 
low altitudes is prevented by both the gas pressure gradient and 
the toroidal component of the magnetic field. This is possible 
because the plasma /3 (J3 - p/(B^/En)) at the surface is rela- 
tively high, varying from 0.16 at the pole to 4.7 at the equator 
Otherwise, in a low /3 corona, the field lines would bend towards 
regions of lower magnetic pressure so that different magnetic 
surface distributions would generate similar coronal topologies, 
as discussed in Sect. |2] Beyond the critical points, the Lorentz 
force that tends to collimate the field lines towards the poles (IX) 
is balanced by the gas pressure gradient (VII), so that the field 
lines remain radial in the poloidal plane. For a significantly more 
magnetic and more rapidly rotating star than the one considered 
here, there is no force capable of balancing the Lorentz force 
and the model cannot generate physically acceptable solutions. 
In other words, the assumptions of no meridional components of 
magnetic field and velocity and of a full open magnetic field is 
incompatible with a magnetically dominated low corona. 

The (^-component of the equation of motion expresses that 
the change of angular momentum of the wind is equal to the 
magnetic torque. As a consequence, the azimuthal velocity is 
close to co-rotation near the surface due to the strong magnetic 
torque, while far from the surface this torque becomes less ef- 
fective and the conservation of angular momentum implies that 
the azimuthal velocity must decrease proportionally to l/R, as is 
illustrated in Fig.|6] 

Figure |2]represents the gas pressure and temperature at two 
different latitudes as a function of radial distance. The pressure 
close to the surface increases towards the equator but has the 
opposite behaviour further out, i.e., as pointed out earlier, the 
gas pressure gradient changes sign close to the Alfven point. In 
this model, the coronal temperature increases towards the pole 
as a consequence of the latitudinal behaviour of the gas pressure 
and density. Having in mind that the temperature is not imposed 
but obtained a posteriori from the perfect gas law, its values of 
T 10^ - lO^K are in relatively good agreement with what we 
would expect for an active solar-type star. It is usually assumed 
in numerical models of solar and stellar winds that the corona 
is isothermal, or, more generally, it follows a polytropic law. In 
our approach it is not possible to impose this behaviour and it is 
also not desirable, as stellar coronae are not well described by 
such simple laws. Additionally, it is not expected that coronal 
regions of large and small magnetic field concentration have the 
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Fig. 4. Absolute values of the different terms of the radial component of the force balance equation as a function of radial distance. 
Positive terms are labelled by (+) and negative terms by (-). In panel a) at 45° of latitude. In panel b) at the equator. 



same temperature. Instead, the temperature profile allows us to 
deduce, from a consistent solution of the energy equation, the 
regions where energy is deposited or removed 

where F is the ratio of specific heats and cr represents the net 
effect of all sources and sinks of energy per unit of mass. In 
the particular case of the solution presented here, the temper- 
ature varies with latitude and with radius. From Figs. |7] and [8] 
we can infer that the temperature and heating rate are higher at 
high latitudes than at low latitudes and that there must be a heat- 
ing mechanism that deposits energy at large radii, beyond the 
sub-Alfvenic region. Although this model has no specific heat- 
ing mechanism included, we note that Alfven waves dissipation 
is a viable mechanism for energy deposition far from the sur- 
face, generating tem perature profile s somewhat resembling the 
ones presented here jCranmeii l2005l) . 

The exploration of the parameter space reveals some impor- 
tant results concerning the overall behaviour of these wind solu- 
tions. In general, an arbitrary choice of parameters of the model 
will not generate wind solutions. In particular for the anisotropy 
parameters, /i and 6 have a crucial influence in the properties of 
the solutions. On one hand, accelerating wind solutions require 
the density to be higher at the equator than the pole, i.e., positive 
values of 6. On the other hand, there are meaningful solutions 
only for fi < 0. If jj takes positive values, the critical solution 
shows a different fast magnetosonic critical point - a spiral type 
point, preventing any possibility of building wind-type solutions. 
This is shown in Fig. |9] where a comparison between a typical 
wind solution for a negative value of fi and a terminated solution 
for a positive value of fj. is presented. 

The effect of the different parameters on the wind solution 
can only be understood by taking into account force balance 
along the lines (radial and azimuthal components of the mo- 
mentum equation, Eqs.[T0land[T2b but also from the equilibrium 
across the lines (0-component of the momentum equation, Eq. 
fTTI ). As an example, if we change the parameters in order to de- 
scribe an even more magnetic and more rapidly-rotating star, we 
would not be able to obtain a solution. Physically, this results 
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Fig. 5. Absolute values of the different terms of the f^-component of 
the force balance equation as a function of radial distance at the pole. 
Positive terms are labeled by (+), negative terms are labeled by (-) and 
terms that change sign by (+) together with a cross along the part where 
they are negative. 



from lack of force balance across the field. The low gas pressure 
can neither balance the tendency of field lines to bend towards 
low latitudes at small distances from the surface, nor the ten- 
dency of field lines to collimate at large distances as required by 
the assumptions of the model. 

A complete study of the behaviour of the wind solution with 
the variation of the parameters is extremely complex because 
the influence each parameter has on the critical solution is of- 
ten dependent on the values taken by the other parameters. This 
particular study is beyond the scope of the present work. 

3.3. Variation of angular momentum loss with the distribution 
of the surface field 

In order to compare the angular momentum loss rate for dif- 
ferent distributions of the magnetic field at the surface of the 
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Fig. 6. Azimuthal velocity at the equator as a function of radial dis- 
tance. Also represented are the profiles corresponding to co-rotation 
(dot-dashed Une) and angular momentum conservation (dashed Une) 
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Fig. 7. Radial profiles of pressure and temperature at 45° of latitude 
and pole. 




Fig. 8. Radial profile of the heating rate per unit of mass at 45° of lati- 
tude and pole. 




Fig. 9. Comparison between solutions for positive and negative values 
of the parameter fx. The solid bold line represents the solution typical 
topology for wind type critical solutions (// = -0.05). The narrow solid 
line represents the topology for solutions with positive values of [i (= 
0.005), which are not wind solutions. The * markers stand for star-type 
points (the Alfven points), the x marker represents the x-type point and 
the s marker points out the spiral-type critical point. 



star, we must prescribe a fixed value of total magnetic flux, Fq. 
Therefore, the surface magnetic field strength at the pole must 
be defined accordingly: 



(14) 



Anrl ^''^ sin 6 yjl + ^ism^' OdO 
The mass loss rate per unit sofid angle at co-latitude 6 is 

- m{e) = pVyT^ = poVorl ^{1 + usin^^ 6X1 + 6sm^' 6). (15) 

The total angular momentum loss rate is equivalent to that car- 
ried by the gas kept in co-rotation with the star out to the Alfven 
surface, which can be expressed as 



- j{e) = An 



1 '■ 



D.(e)p,v,Rts^ Ode. 



(16) 



Or equivalently, in terms of the model parameters (Lima et al., 
2001) as 



j(0) = Ar. 



sin"+2 o^\+n sin^' 0d0. 



(17) 



In this model, p, Vr and Q are functions of 5 and yield an an- 
gular momentum loss rate that is not directly dependent on this 
anisotropy parameter. By changing the anisotropy parameters ju 
and e the Alfven radius and Alfven velocity change. 

In order to evaluate the polar concentration of a given mag- 
netic field surface distribution, we first compare the magnetic 
flux at low co-latitudes (from 30° towards 0°) with the overall 
total flux. Thus, 



-^l+psin^' Ode 



(18) 



Then, we define the parameter t that evaluates the magnetic 
field concentration towards high latitudes by comparing it with 
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Fig. 10. Total angular momentum loss variation as a function of polar 
field concentration, for difi'erent sets of /j (dashed lines) and 6 (solid 
lines). 

a uniform magnetic field distribution (the split-monopole case 
with e = 1 and jU = 0) 

T= =3^^^^ ^ (19) 

'^monopole ^1 +jjsm^'0de 

If T > 1 the field is concentrated towards the pole while if 
T < 1 the field is concentrated towards the equator 

As a result of the way the model is constructed, a change 
in the parameters /i or e will change the total amount of mag- 
netic flux, the location of the Alfven radius and also the rotation 
period of the star (see Eqs. [14] and |9]l. Therefore, an iterative 
procedure on the parameters A and M4 is performed so that the 
different solutions under comparison have the same total amount 
of magnetic flux and rotation period within an eiTor of 5%. 

By evaluating the angular momentum loss rate for different 
wind solutions attained for several sets of fi and e, we are able 
to obtain a relation between the magnetic polar field concen- 
tration and the total angular momentum loss. This is shown in 
Fig-Hni where solid lines represent how - J varies with t due to 
the variation of fj., with constant e. Physically, this coiTesponds 
to determining how the variation of the magnetic field distribu- 
tion with a fixed rotation profile affects the angular momentum 
loss rate. Also represented by dashed lines is the case of how 
-J varies with t due to the variation of e, with constant fi. In 
this case, both the magnetic field and rotation latitudinal profiles 
vary. All solutions considered here are lied as SMR, revealing 
that centrifugal forces are not predominant for the wind accel- 
eration mechanisms. The dynamics of the wind is, as previously 
shown, independent of the forces generated by the rotation of the 
central body. This is true even for the particular case where the 
period is of 4 days since other contributions, such as the pressure 
gradient, are more important. 

Our purpose is to study young and rapidly rotating stars, 
therefore we must consider small values of e and small values of 
6, as this kind of object shows near rigid body surface rotation 
with a slightly faster equator (e.g. Collier Cameron & Donati, 
2002). Also, it is only adequate to compare stars with differ- 
ent surface magnetic field distribution but with the same surface 



rotation latitudinal profile. Thus, we restrict the physical appli- 
cation of this model to the cases of constant e and variable fi. If 
we analyze the solid lines in Fig. (TO] we see that the total an- 
gular momentum loss is an increasing function of r. Therefore, 
the higher the magnetic field polar concentration, the more effi- 
cient magnetic braking is. Referring to Eq.[T6l and noticing that, 
by continuity, p^VtR* = poVo/Jy/jJ, changes in total angular mo- 
mentum loss must arise from changes in the mass loss rate or 
changes in the Alfven radius. But as the mass loss rate decreases 
with increasing r and \fi\ (cf. eqfTSll, we can conclude that the 
increase in total angular momentum loss is due to an increase 
in the Alfven radius that more than compensates the decrease in 
mass loss. 

As |yu| increases, Bq increases and decreases, i.e., at the 
pole the star becomes more magnetic as the field concentrates 
there. As V,RI oc (M^)-^ we see that either the Alfven velocity, 
or the Alfven radius or both must increase. From the different so- 
lutions present in Fig.[TO]as well as an analysis of the differential 
equation yielding Vr{R) (or Y(R)}, we find that the wind acceler- 
ation decreases as fi increases so that y» decreases and the Alven 
radius must increase. Physically, this increase in the Alfven ra- 
dius is due to a less efficient acceleration mechanism. The anal- 
ysis and clear understanding of this result is complicated by the 
fact that as yU changes there is a direct effect on the radial mag- 
netic field but also on the radial velocity latitudinal profile and, 
indirectly, the wind pressure and temperature. Therefore, the in- 
crease in angular momentum cannot be solely attributed to the 
increase of magnetic field concentration as this cannot be sepa- 
rated from changes in other physical quantities. This result is at 
odds with the result of Sect. |2] where no significant difference in 
angular momentum loss as a function of field concentration was 
found. This discrepancy may be explained by the fact that the 
two models address totally different scenarios characterized by 
different plasma /3 regimes. 

4. Discussion 

Empirical or qualitative models of magnetized stellar winds have 
been very successful in establishing the basic physical prin- 
ciples that govern the rate of angular mo mentum l oss a ssoci- 
ate d with a wind. In pa rticular, the work of iMestell (11968 ) and 
Mes tel & Spruit! (Il987l) investigates how the magnetic braking 
efficiency varies with stellar rotation rate. It shows that it is pos- 
sible to have saturation in the angular momentum loss rate with 
rotation without dynamo saturation as a result of two competing 
effects. On the one hand, the magnetic torque increases with ro- 
tation rate due to the increase of the field strength. On the other 
hand, the fraction of the stellar surface with open magnetic field 
lines contributing to braking decreases with rotation rate. This 
model has also been applied to establish that a rearrangement 
of the surface field from a low order to a high order multipole 
at high rotation rates implies a decrease in braking efficiency 
(lTaam&Spruillll989l) . 

There are, however, limitations to how much information can 
be obtained from these qualitative models. In this work, we sim- 
ply address the question of whether a correct estimate of the in- 
fluence of the magnetic surface field concentration towards the 
poles can be obtained using these empirical models or whether 
such a goal requires a computation of the field that takes into 
account force balance in all directions. To this end, in Sect. |2] 
we determine the coronal magnetic field due to a surface flux 
distribution concentrated around the poles and show that this 
distribution is not maintained further out. Instead, the magnetic 
field rapidly expands to regions of low magnetic pressure and ap- 
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proaches the simple split-monopole field. Flux tubes at different 
latitudes have very different expansion rates, but we find that this 
has no significant effect on the angular momentum transported 
outwards. 

In order to investigate the properties of stellar winds for dif- 
ferent surface flux distributions we apply, in Sect. [3] the analyt- 
ical wind model of Lima et al. (2001). To our knowledge, this 
is the first application of an analytical quantitative 2D model to 
the problem of magnetic braking in young solar- type stars. This 
model has very appealing characteristics has it describes the out- 
flow of a rotating star for which the latitudinal magnetic flux 
distribution is the same at all radial distances. It also has some 
limitations that are typical of this kind of model. It is constructed 
under the assumption of separation of variables and there is no 
energy equation, with the temperature being determined from 
the perfect gas law. Also, although the model allows different 
surface flux distribution, its range of variation is somewhat lim- 
ited. Its most restrictive properties, however, are that the Alfven 
surface is spherical, contr ary to what is expecte d for fast rotators 
and has been proposed bv lSolanki et alj (Il997l) . and the absence 
of meridional components of the velocity and magnetic field. 
Therefore, the conclusions obtained from using this model can- 
not be a priori considered as general features. Yet, some of the 
results are extremely relevant as they convey fundamental physi- 
cal principles. We find that the surface concentration of magnetic 
flux at high latitudes can only occur in the corona for a relatively 
high plasma beta, so that the gas pressure gradient can oppose 
the Lorentz force that tends to bend the field lines towards the 
equator. This work stresses the importance of the latitudinal pro- 
files of the physical quantities involved in the description of the 
hydro-magnetic wind. As shown in section 13.21 the radial fea- 
tures of the wind, such as the acceleration mechanisms, are not 
only dependent on the force balance along the lines but also on 
the force balance across the lines. Furthermore, a given radial 
field latitudinal profile implies, due to force balance and under 
the assumptions of the model, a certain latitudinal profile of the 
density and velocity that are as important as the magnetic field 
in determining the angular momentum loss rate. The most im- 
portant prediction from this model is that a higher polar field 
concentration leads to larger braking rates than a smoother field 
distribution. However, because of the inherent coupling of this 
result with some of the assumptions of the model, this cannot 
be viewed as general. In spite of this, we have presented an ex- 
ample of an equilibrium model, with unknown stability proper- 
ties, for which magnetic braking increases as the field concen- 
tration to the pole increases. Remarkably, this model has the at- 
tractive pr operty of ha ving a magnetic topology ide ntical to the 
models of [Holzwarflj (2005), iHolzwarfli & Jardind (2006) and 
ISolanki et all (1997) but it complies with force balance in all 
directions, which generates antagonistic results. 

Both models presented in this work neglect the poleward col- 
limation of the field lines characteristic of FMR. An interesting 
question is whether this poleward collimation increases or de- 
creases the angular momentum loss of the star, but we are not 
aware of any work specifically addressing this issue. At high lat- 
itudes increases due to the effect of collimation but the op- 
posite happens near the equator (e.g. lSakuraii 1985), so that the 
end result depends on which effect is dominant. In Sect. |2] we 
have determined that at 4 - 8 stellar radii from the stellar sur- 
face the radial field becomes uniform in latitude, independently 
of its surface distribution, and that in general this occurs very 
much inside the Alfven surface. Therefore, we expect that field 
collimation would occur identically for both uniform and highly 



concentrated surface field distributions with no significant dif- 
ferences in angular momentum loss. 

Based on the results presented in Sect. |2]&|3] we argue that 
the concentration of magnetic flux at high latitudes does not di- 
rectly contribute to limit the braking efficiency of the wind. By 
comparing the results obtained with the two different models one 
could in principle be able to infer how crucial the limitations 
introduced are in the analytical treatment of the 2D model to 
the results obtained. The extent to which this comparison can 
be made is limited by the very different assumption about the 
plasma [} in the two cases. We hypothesize that the increase in 
angular momentum loss with magnetic field concentration ob- 
served in Sect. [3] results from changes in the gas temperature, 
density and pressure that are generated to maintain force bal- 
ance across the field in the absence of meridional components. 
One can speculate about what to expect if some of the most strin- 
gent assumptions of our 2D model are relaxed. We suggest that 
allowing for meridional components, non spherical Alfven sur- 
face and a more realistic plasma yS would lead to closed magnetic 
field lines near the equator Furthermore, a surface field concen- 
trated towards the pole would have a weaker field at low latitudes 
and consequently have more open magnetic flux than a smooth 
surface magnetic field. This would again imply larger braking 
rates for a surface field concentrated towards the pole than for 
a smooth field with equal amounts of magnetic flux, if the wind 
properties are similar in the two cases. However, the existence of 
large amounts of magnetic flux at high latitudes can still lead to a 
reduced angular momentum loss rate if it creates a complex field 
topology akin to multi-order magnetic fields (cf . iTaam & Spruitl 
1989). 

The work presented here can be extended by further analyti- 
cal modelling or numerical simulations. However, many aspects 
are still unknown and this constrains how much one can pre- 
dict about magnetic braking. It is still not known whether polar 
spots ar e largely unipolar , as assumed here (for op posing views 
see e.g. lSchriiver&Titlel 2001 and Mclvor et al.L 2003). In ad- 
dition, we cannot rule out the possibility that stellar winds are 
accelerated by dissipation of magnetic waves as well as thermal 
and centrifugal forces. The nature of the slow solar wind is still 
largely unknown and it remains an open question whether it is 
admissible to ignore the contribution of the slow wind to the an- 
gular momentum evolution of active late-type stars. 



5. Conclusions 

In the present work we investigate how the magnetic surface field 
distribution affects the coronal magnetic field and the rate of an- 
gular momentum removal by the stellar wind. 

There are three important results from our work: First, we 
have shown that very different surface flux distributions yield 
similar coronal fields in a low yS coronal plasma as well as similar 
wind braking rates. Second, we have demonstrated that the radial 
features of the wind, such as the acceleration mechanisms and 
the gas pressure distribution, are not only dependent on the force 
balance along the field lines, but also on the force balance across 
the field lines. Finally, in the wind model of (Lima et al., 200 1]) 
a higher polar field concentration leads to larger braking rates 
than a smoother field distribution. However, we cannot rule out 
the possibility that this is a result of the assumptions of the model 
and so it may not be regarded as a general feature. This model 
also demonstrates that the rate of braking is dependent on the 
latitudinal behaviour of several physical quantities and not only 
of the surface radial field. 
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We conclude that the concentration of magnetic flux at high 
latitudes is unlikely to directly constrain the braking efliciency 
of the wind. It may, however, have a decisive importance in de- 
termining the amount of open magnetic flux that contributes to 
wind magnetic braking. 
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Appendix A: Stream functions 

Here we give the explicit forms of Z„ for n - 5 and n -1: 
15 
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Appendix B: Wind model 

Here we present how the wind solution is obtained. As we fol- 
low very closely the approach of Sakurai ( 1985), only a brief de- 
scription is given. We assume azimuthal symmetry and by com- 
bining the equations of mass and magnetic flux conservation, the 



frozen-in condition, the polytropic law and the momentum con- 
servation equation in the poloidal and azimuthal directions we 
arrive at a Bernoulli integral of the equation of motion 



2 



r p 



y-lp 



2 
GM 
r 



^ p. 



(B.l) 



where Vp is the poloidal velocity, Q the stellar angular velocity, 
M the mass of the star, G the gravitational constant, E an inte- 
gration constant and the subscript * stands for the Alfven point. 
Upon applying the law of mass conservation, pVpS - p^)Vp^^S^), 
with S representing the area of the flux tube and the subscript 
the values at the stellar surface, we obtain H - H(r,p). At the 
fast and slow critical points one has the regularity conditions 



dH _dH _^ 
dp dr 

together with 

H(r„p,)^H(rf,pf)^E. 



(B.2) 



(B.3) 



Applying mass and flux conservation allows us to write the den- 
sity at the Alfven point as 



p. =Po 



(B.4) 



The system of six simultaneous algebraic equations ( IB.2IIB.3l l is 
solved for the unknowns r,, r,, p,, rj, pf and vo. The integration 
constant E is not an additional unknown as it can be w r itten i n 
terms of Vo and r». We note that, contrary to ISakurail ( Il985h . 
we solve our equations in dimensionless form with respect to 
surface values and not with respect to the Alfven radius as, in 
general, S cc does not hold. 

References 

Altschuler M.D., Newkirk, Jr. G.: 1969, Sol. Phys. 9, 131 

Balogh A., Smith E.J., Tsurutani B.T., et al.: 1995, Science 268, 1007 

Barnes J.R., Collier Cameron A., Unruh Y.C., Donati J.F, Hussain G.A.J.: 1998, 

MNRAS 299, 904 
Barnes S., Sofia S.: 1996, ApJ 462, 746 
Belcher J.W., MacGregor K.B.: 1976, ApJ 210, 498 
Buzasi D.L.: 1997, ApJ 484, 855 

ColUer Cameron A., Donati J.F: 2002, MNRAS 329, L23 
Collier Cameron A., Robinson R.D.: 1989, MNRAS 236, 57 
Cranmer S.R.: 2005, In: ESA SP-592: Solar Wind 11/SOHO 16. Connecting Sun 
and Heliosphere 

Donati J.F, Collier Cameron A., Hussain G.A.J., Semel M.: 1999, MNRAS 302, 
437 

Heyvaerts J., Norman C.A.: 1989, ApJ 347, 1055 

Holzwarth V.: 2005, Ab-A 440, 41 1 

Holzwarth V.. Jardine M.: 2005, AM 444, 661 

Holzwarth V., Jardine M.: 2006, ArXiv Astrophysics e-prints 

Hussain G.A.J., Jardine M., Collier Cameron A.: 2001, MNRAS 322, 681 

Keppens R., Goedbloed J.R: 1999, A&rA 343, 251 

Lepeltier T., Aly J.J.: 1996, AM 306, 645 

Li J.K., Wu K.W., Wickramasinghe D.T.: 1994. MNRAS 270. 769 
Lima J.J.G., Priest E.R., Tsinganos K.: 2001, A6-A 371, 240 
Low B.C.: 1986, ApJ 310, 953 
MacGregor K.B., Brenner M.: 1991, Ap/ 376, 204 

Mclvor T, Jardine M., Cameron A.C., Wood K., Donati J.F: 2003, MNRAS 345, 
601 

Mestel L.: 1968, MNRAS 138, 359 

Mestel L., Spruit H.C.: 1987, MNRAS 226, 57 

Parker E.N.: 1979, Cosmical Magnetic Fields. The International Series of 
Monographs on Physics, Oxford University Press, New York 



A. Aibeo et al.: Polar spots and magnetic braking 



11 



Riley P., Linker J.A., Mikic Z., et al.: 2006, ApJ 653, 1510 

SakuraiT.: 1985, A&A 152, 121 

Sauty C, Tsinganos K.: 1994, AfrA 287, 893 

Schatzman E.: 1962, Anjia/e.s- d'Astwphysique 25, 18 

Schmitt J.H.M.M., Favata E: 1999, Nature 401, 44 

Schrijver CJ., Title A.M.: 2001, ApJ 551, 1099 

Schuessler M., Solanki S.K.: 1992, A&rA 264, L13 

Smith E.J., Balogh A.: 1995, Geophys. Res. Lett. 22, 3317 

Solanki S.K., Motamen S., Keppens R.: 1997, A&A 325, 1039 

Strassmeier K.G.: 2002, Astronomische Nachrichten 323, 309 

Strassmeier K.G.: 2005, Astronomische Nachrichten 326, 269 

Suess S.T., Richter A.K., Winge C.R., Nemey S.E: 1977, ApJ 217, 296 

Suess S.T., Smith E.J.: 1996, Geophys. Res. Lett. 23, 3267 

Taam R.E., Spruit H.C.: 1989, ApJ 345, 972 

Tsinganos K., Sauty C, Surlantzis G., Tnissoni E., Contopoulos J.: 1996, 

MNRAS 283, 811 
Tsinganos K., Trussoni E.: 1991, A&A 249, 156 
Vogt S.S., Penrod G.D.: 1983, PASP 95, 565 
Weber E.J., Davis L.J.: 1967, ApJ 148, 217 



